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ABSTRACT 

The puzzles of black hole evaporation can be studied in the simplified context of 
1+1 dimensional gravity. The semi-classical equations of Callan, Giddings, Harvey 
and Strominger provide a consistent description of the evaporation process which 
we describe in detail. We consider the possibility that black hole evolution leads to 
massive stable remnants. We show that such zero temperature remnant solutions 
exist but we also prove that a decaying black hole cannot evolve into one of them. 
Finally we consider the issue of loss of quantum information behind the global 
event horizon which develops in these geometries. An analogy with a well known 
solvable system shows that there may be less to information than meets the eye. 
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1. Introduction 



Black hole evaporation provides us with an interesting set of puzzles where 
quantum mechanics and general relativity clash [1,2,3]. The central paradox is 
that matter in a quantum mechanically pure state can undergo gravitational col- 
lapse to form a black hole, which then evaporates, and in this process the initial 
pure state appears to evolve into a mixed quantum state, describing the outgoing 
Hawking radiation. This argument is based on a calculation which treats matter 
quantum mechanically in a background geometry which evolves according to Ein- 
stein's classical field equations. It seems prudent to include the perhaps subtle 
effects of back-reaction on the geometry, due to the emission of Hawking radiation, 
before we start tinkering with quantum mechanics. Hopefully these issues can be 
resolved without direct reference to Planck scale physics (after all the horizon of 
a massive black hole forms in a region of weak gravitational coupling) allowing us 
to explore interesting features of quantum gravity without having to grapple with 
a theory of everything. 

It is often illuminating to address difficult physical problems in the simplified 
context of 1 + 1 dimensions. A particularly interesting toy model for black hole 
physics was suggested by Callan, Giddings, Harvey and Strominger (CGHS) in 
reference [4], where they coupled two-dimensional gravity to a dilaton and con- 
formal matter fields. CGHS described a set of exact classical solutions, including 
black holes formed by infalling shock waves, and they observed that Hawking ra- 
diation can be described by adding correction terms to the classical equations 
of motion. These extra terms derive from the conformal anomaly of the matter 
fields [5] and are local in conformal gauge, with the usual non-locality appearing 
only in the form of boundary conditions required to satisfy constraints. CGHS also 
pointed out that once these quantum corrections have been included one can study 
the back-reaction on the geometry, due to the Hawking radiation, using classical 
methods in the effective theory. They argued that the quantum corrections would 
prevent the usual black hole singularity from forming in gravitational collapse and 
therefore quantum information would not be lost. Subsequently it was shown that 
in this theory gravitational collapse always develops a singularity [6,7]. Thus the 
paradox has not been resolved. 

In this paper we study the CGHS model further. Background material can be 
found in references [4], [6] and [7]. We begin in section 2 by briefly reviewing the 
model and describing some exact classical solutions. We then turn to the quantum 
theory in section 3. Following CGHS we study Hawking radiation and its back 
reaction on black holes using effective equations of motion, with correction terms 
induced by the conformal anomaly of the matter fields. Using both analytic and 
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numerical methods we are able to present a fairly detailed picture of the formation 
and quantum mechanical evaporation of 1 + 1 dimensional black holes. 

An important issue to settle is the nature of the final state that a 1+1 dimen- 
sional black hole will settle down to. One can think of a number of possibilities, all 
of which are problematic to some degree. The first, suggested by Hawking [2], is 
that the black hole completely evaporates leaving no remnant behind but only in- 
coherent radiation, and quantum coherence is lost in the process. This is a radical 
proposal and more conservative alternatives should be exhausted before we modify 
quantum mechanics. 

Another possibility is that the Hawking process leaves behind a stable remnant 
whose precise quantum state reflects the initial configuration of the system [8]. 
Evidently this requires a virtually infinite number of distinct quantum states of 
black hole cadavers. There are two very different cases to consider. On the one 
hand the mass and size of the remnants would be of order the Planck scale, in which 
case their density of states must be infinite. If these light remnant states couple 
to gravity in any conventional way, they would be copiously produced in pairs 
whenever there is available a Planck scale amount of energy. The energy would 
not have to be very concentrated because even if the probability for producing a 
particular pair of these Planck scale objects is tiny the infinite density of states 
makes the total rate diverge. For example, the sun would have been short lived 
indeed if such channels had been available. It has recently been suggested that the 
black hole remnants cannot be produced in pairs [6] but this in itself is a radical 
departure from the principles of relativistic quantum mechanics. 

Alternatively the Hawking process might terminate while the black hole still 
has macroscopic mass. That could only happen if the back-reaction on the geom- 
etry develops a singularity on the apparent horizon itself. This may sound out- 
landish but in section 3 we will show that the curvature on the apparent horizon 
does in fact increase as the black hole evaporates and in section 4 we exhibit static, 
non-radiating solutions of the CGHS-equations which are candidates for such rem- 
nants. Massive remnants can carry a lot of information and their density of states 
would not have to diverge in order to preserve quantum coherence, so this possi- 
bility seems worth exploring. However, we are able to prove that it is not realized 
in 1+1 dimensional gravitational collapse, as described by the CGHS-equations. 

Perhaps the most conservative option is that despite its thermal appearance 
Hawking radiation really carries off information about the quantum state of the 
black hole [3] . This would reconcile the evolution of black holes with conventional 
quantum mechanics. In section 5 of this paper we will discuss to what extent this 
type of questions can be addressed in this 1+1 dimensional model, but unfortu- 
nately we cannot offer any definitive conclusions at present. 
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2. The Model 



The system we will consider is two-dimensional gravity coupled to a dilaton <j) 
and a set of scalar fields f\. The classical dynamics is governed by the action 

1 r 1 N 

S = — J d 2 x V^[e- 24> (R + 4 (V0) 2 + 4A 2 ) - - E( V ^)'] » ^ 

which resembles the target space action of c—1 non-critical string theory. It can be 
viewed as an effective action for radial modes of four-dimensional extreme dilaton 
black holes [4,9,10] but in this paper we will study it on its own merits as a two- 
dimensional model of gravity coupled to matter. 

The matter fields fi are minimally coupled to the metric and do not directly 
interact with the dilaton. We will assume there is a large but finite number of 
them. The reason for this is that the one-loop quantum corrections to the action 
are dominated by the conformal anomaly of the matter fields if N is large and 
the effective action takes a particularly simple form. The arbitrary parameter A 
defines a mass scale in this theory and we will refer to a black holes with M > NX 
as massive. 

The classical equations of motion obtained from the action (2.1) are 

v 2 /, = o, 

R + 4A 2 + 4V 2 - 4(V(/>) 2 = , 



1 N 

2e~^ [V M V,0 + ^,((V0) 2 - V 2 - A 2 )] =- X)(V M /iV„/i - ^(V/ 4 ) 2 ) 



i=l 

(2.2) 



It is convenient to work in conformal gauge 



9+- =9-+ = ~ \e 2p 
9— = 9++ = , 



(2.3) 



with light-cone coordinates x — x ± x . The fields in the theory are then fi, 
and p, and the classical equations of motion become, 

d+d-f t =0, 



A_ 2 
2 

A 2 



2d+8-4> - 2d + 4>d-4> - -e 2 P = d+d-p , (2 . 4) 



d+d-<f> - 2d+(f>d-(f) e 2p = , 

respectively. In addition, one must impose as constraints the equations of motion 



4 



corresponding to the metric components, which are set to zero in this gauge: 



(2-5) 



2 

i=i 



A number of exact solutions of these equations can easily be found [4]. The 
simplest one is the vacuum solution 



fi = 0, 

™ 9 9 ^ (2-6) 

e -20 = e -2 P = _ A 2 x+x - 



This is just the 'linear dilaton' background of non-critical string theory. A change 
of variables, x ± = ±e ±u , makes the metric flat and the dilaton field linear in 
u + —u~. A 'static' black hole is described by 

/,: =0 , 
A 

Finally one can describe the metric and dilaton fields due to an infalling shell of 
matter by patching together the vacuum solution and a black hole solution across 
some light-like line x + = Xq (see figure 1), 

{— \ 2 x + x~ if x + < xt\ 

(2%) 
{ x +- x +)-X 2 x+x- i£x+>x+. [ } 

Note that (j) and p are continuous across the matching line. Inserting (p and p from 
(2.8) into the ++ constraint equation in (2.5) reveals that this solution is indeed the 
response to an incoming shock wave with energy-momentum ^ 2^2iLi d+fid+fi = 
j^d{ x+ — x o)- The geometry is singular on the curve where e~ 2p = and the 
line x~ = — rr~T is a global event horizon in this classical solution. 

The factor of e _2< ^ in front of the gravity terms in the action (2.1) makes the 
strength of gravitational quantum corrections depend on the magnitude of the 
dilaton field, g ~ e^. Both the linear dilaton vacuum and the classical black hole 
solutions have a position dependent coupling strength with weak coupling being 
asymptotically far away. In fact, the strength of quantum corrections provides a 
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coordinate invariant notion of position in this one- dimensional world. Recall that in 
the four-dimensional Schwarzschild solution lines of constant radial distance (and 
fixed angle) go from being time-like outside the horizon to being space-like inside 
the black hole. The analogous statement for our 1+1 dimensional black holes is 
that (V0) 2 changes sign at the horizon [7]. This provides us with a sensible local 
definition of an apparent horizon which will be useful in what follows. 



3. Quantum Corrections 

The most important result of the original paper of CGHS [4] is that the back- 
reaction on the geometry due to Hawking radiation can be described by adding 
appropriate correction terms to the classical equations of motion. When there is a 
large number of matter fields in the theory there is a class of quantum corrections, 
associated with closed matter loops, which become important when g 2 N ~ 1, long 
before other gravitational corrections kick in. At the one-loop level the matter 
conformal anomaly is the only term of this class and it will be responsible for the 
quantum back-reaction we will explore later on. When a black hole is formed the 
coupling at the horizon is c/ 2 ~ so our approximations remain valid so long as 
the mass remains large compared to NX. 

At one-loop order the anomaly contributes the well known non-local term 

1 ' d 2 xy/=gR^R (3.1) 



96vr / y V 2 



to the effective action [11]. In the conformal gauge (2.3) this term appears local, 
— TTk I d 2 xd+pd-p and the one-loop corrected equations of motion, which we refer 
to as the CGHS-equations, can be arranged as follows, 

d+d-fi=0, 

d + d-<P = (l~e^)d+d-p, ( 3. 2) 
2 (1-^) d+d-tp = (1-^e 2 *) (4«9 + 0d_0 + A 2 e 2 ") . 



★ Following CGHS we assume that the Liouville cosmological constant term can be fine tuned 
to zero. 
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The constraint equations also receive corrections, 



1 N N 
^(2d 2 - 4d +P d+<P) = - d+fid+fi - - (8+pd+p - d 2 +P + t+(x+)) 



1 2 ' •" 12 

i=i 

N 

-20 



e -(2«9 2 - 4d_pd_0) = 1 ^d- fid- ft - ^(d_pcLp - «9 2 p + t-(x~)) . 

i=i 

(3.3) 

The functions t±(x ± ) are needed to satisfy asymptotic physical boundary condi- 
tions [4], so there remains a degree of non- locality in the conformal gauge. As 
expected the quantum corrections in these equations appear with factors of Ne 2 ^ 
compared to the classical terms. 

For some purposes it can be useful to have a description of this system outside 
conformal gauge. In a more general gauge the anomaly induced term (3.1) can 
be made local by adding a new degree of freedom which is coupled to the scalar 
curvature, 

i r 1 / at 

(3.4) 



The non-local form (3.1) is recovered when Z is integrated out of the path integral. 
The energy momentum tensor of the Z field is 

V = l - [V^ZVvZ - ^(VZ) 2 ] + yj[V M V,Z - g^V 2 Z] . (3.5) 

The conformal anomaly of the matter fields appears as a non-vanishing trace of 
the classical energy-momentum tensor of Z. This is analogous to what happens in 
bosonization of two-dimensional spinor electrodynamics where the axial anomaly 
of the fermions can be expressed classically in terms of the bosonizing field. We 
will develop this analogy further in section 5 when we discuss the fate of quantum 
information in this theory. 

The classical equations (2.4) can be made linear by a combination of gauge 
choices and field redefinitions, which makes them relatively easy to solve. Un- 
fortunately these tricks no longer work when the one-loop corrections have been 
added, and we now have coupled, non-linear, partial differential equations to deal 
with. It is easily checked that the linear dilaton vacuum (2.6) remains an exact 
solution. The geometry due to an incoming shock wave can no longer be found in 
closed form but is described by some non-trivial solution of the CGHS-equations 
which is matched continuously onto the vacuum across x + = Xq. It approaches a 
classical black hole solution asymptotically far away, where the coupling is weak. 
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In the remainder of this section we will use the equations of motion and con- 
straints to explore qualitative features of such solutions, reinforcing and extending 
the observations of reference [7]. We will also present numerical calculations of 
two-dimensional black hole evolution which support our analytic results. 

The first thing to note is that a curvature singularity forms on the infall trajec- 
tory at e 2 * = This was shown in [6] and [7] and we will not repeat the argument 
here. The singularity extends into a curve in the x + > Xq region through which 
physical solutions cannot be meaningfully extended. Our approximations break 
down as the singularity is approached and it may or may not persist in the full 
quantum theory, but its presence at this level means that care must be taken in 
developing any systematic approximations. Perhaps one should view the y^e 2 ^ = 1 
line as a natural boundary to spacetime in the quantum theory. We will at any 
rate focus on the evolution of massive black holes where the CGHS-equations are 
reliable so long as the mass remains large compared to iVA. Our results are for 
the most part independent of how the singularity issue is ultimately resolved. It is 
worth pointing out that removing the singularity by modifying the short distance 
physics will not automatically cure problems of information loss. The evapora- 
tion process would still be governed by the semi-classical equations for most of the 
black hole lifetime and there will be a limited amount of energy left to carry off the 
information contained in the original black hole (we are assuming for the moment 
that it does not get radiated away during the semi-classical evaporation). Either 
it would be contained in a light stable remnant or have to be emitted in very low 
energy radiation over a long time. In either case we would have to contend with a 
virtually infinite density of (meta-) stable states at the Planck scale in the quantum 
theory. 

The singularity is cloaked by an apparent horizon which forms at x~ ~ — , 
for M large, and then begins to recede as the black hole emits energy in Hawking 
radiation [7]. Later on we will show that the apparent horizon continues to recede 
as long as there is outgoing energy flux from the black hole. This intuitively obvious 
result provides a nice consistency check on the picture of black hole evaporation 
presented by the CGHS-equations. 

The picture of black hole evaporation which emerges is qualitatively as follows. 
The curve of singularity at y^e 2 ^ = 1 forms inside the apparent horizon where V</> 
is time-like and lines of constant <fi are space-like. Thus it follows that the line of 
singularity is everywhere space-like unless it were to cross the apparent horizon. 
Such a crossing, however, can not happen since d + (p equals zero on the apparent 
horizon but diverges on the singularity. Once the black hole has formed it emits 
Hawking radiation and the apparent horizon recedes. Since it cannot cross the 
singularity they must both approach light-like lines. We will present analytic and 
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numerical evidence which strongly suggests that they in fact approach the same 
light-like line from opposite sides as in figure 2. The light-like asymptote is a global 
event horizon. 

If all the mass of the black hole is radiated then one expects the geometry 
to approach the vacuum solution. Alternatively the black hole might approach 
a massive stable remnant whose properties would emerge at asymptotically late 
times. To investigate the asymptotic behavior of black hole solutions one can 
consider how coordinate invariants, such as the curvature scalar or (V0) 2 , evolve 
along contours of constant coupling strength. One can, for example, derive a simple 
relation between the rate of change of R along a curve of constant and the flux 
of energy across the curve. For this purpose it is convenient to introduce local 
coordinates u, v in the neighborhood of the contour, such that u parametrizes the 
contour, a; - (it)) = 4>o, and v parametrizes orthogonal curves. A good 

choice is to take v as itself and normalize it in such a way that the Jacobian 
determinant of the coordinate transformation from x + , x~ to it, v equals one. The 
matter energy, including Hawking radiation, in this coordinate system is T uu and 
the outgoing energy flux is T u v . 

The curvature scalar in the conformal gauge (2.3) is given by 

R = 8e- 2p d+d-p, (3.6) 
which can be reexpressed using the equations of motion (3.2) as 

R = ,_Z 2 M e- 2p d + <pd^+\ 2 ) . (3.7) 

1 12 e 

We are interested in the rate of change of curvature along the <j) = 4>o contour, 

R ~du R U=<t>o 

dx+ dx~ (3.8) 

=— — d + R+ — — d-R v > 

ou ou 

=d-<pd+R- d+<pd^R. 
Inserting the expression (3.7) for the curvature gives 

R = i Tk^o [(d-^^-Ad+pd^) - (a + 0) 2 (2^ ^-Ad-pd^)] . (3.9) 

1 12 e ° 

The constraint equations (3.3) reveal that e^ 2 ^(2d\4>—Ad±pd±(j)) is the gravita- 
tional response to the energy-momentum tensor of matter and Hawking radiation 
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so we can write 



-2p+20 o 

k = t-K 2*o T++ - ^ ■ (310) 

1 12 e 

The final step is to observe that the expression in the square brackets is proportional 
to the energy flux across the contour and we obtain the desired result 

A P 2( t>o 
1 12 e 

This relationship can be checked by considering a contour in the classical region 
Ne 2< ^° « 1. The value of the curvature of the classical solution (2.8) where the 
contour meets the infall line is R — A\e 2 ^°M. If the curvature on the contour 
goes to zero at late times then (3.11) implies that the integrated flux of Hawking 
radiation carries off an amount of energy equal to the original mass. In other 
words, the geometry approaches flat space if and only if the black hole radiates 
away all its energy. 

We are now ready to show that the apparent horizon always recedes if there is 
any outgoing flux of energy from the black hole. Let the line of apparent horizon 
be parametrized (x + ). By definition d + (j) vanishes everywhere along 

this curve, 



d 



j2 



dx 



= — d + _ = + (— ) d + d^ . (3.12) 



The constraints (3.3) relate the value of d\<p on the apparent horizon to the out- 
going energy flux at that point. Using this and the second equation of motion in 
(3.2) we can write the slope of the horizon curve in terms of the flux as follows, 

dx+ 2A 2 (l-^)(d_ ( />) 2 u ' 1 ' ) 

The horizon recedes at a rate proportional to the outward flow of energy from the 
black hole (and advances when the black hole absorbs infalling matter). 

On the apparent horizon the expression for the curvature (3.7) reduces to 
R = ^_N_ e 2^ an d it immediately follows that the horizon curvature increases as 

the black hole evaporates and the apparent horizon recedes to larger values of e 2 ^ . 

The apparent horizon of a massive black hole forms in a region of weak coupling, 
e 2 ^ ~ jjj, and a good approximation to ^Lp j s obtained by inserting the classical 
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solution (2.8) on the right hand side in (3.13). To leading order this gives [7] 



dx~ N 
dx+ 48A 2 x+ 



2 , (3-14) 



which corresponds to a rate of energy loss of in agreement with the asymptotic 
flux of Hawking radiation found by CGHS [4] . So long as the remaining black hole 
mass is large compared to NX the rate of energy loss is slow compared to the mass 
and (3.14) (with Xq replaced by x + ) gives the instantaneous slope. Integration 
gives an 'adiabatic' approximation to the horizon curve, 



x 



In the x + — > oo limit this approaches a global horizon at x = x + 48 j[^ + ■ For a 
very large initial mass this should be an excellent estimate because (3.14) is valid 
for most of the black hole lifetime. The approximation breaks down when the 
coupling strength becomes too large. Adiabatic arguments similar to those given 
above lead to a linearly growing coupling strength at the apparent horizon, 

which remains small for a long time if the original mass is large. 

Note that in the coordinate system the total distance along x~ , which the 
apparent horizon recedes, is independent of the black hole mass. This means that 
the global event horizon extends far into the weak coupling region and at first 
sight it appears that some, if not almost all, of the information about the quantum 
state of the system will be lost behind it. However, as we shall soon see, the 
back-reaction is strong enough to convert all the energy of the infalling matter 
to outgoing Hawking radiation by the time the incoming shock wave reaches the 
global event horizon. We will address the question of information loss in section 5 
and at that point it will be useful to have a clear notion of the total energy carried 
across light-like lines of constant x~ , such as the global horizon. 

We will define the total energy on a given x~ line as 

M{x~)= lim ^(l-^)§e- 2 ^. (3.17) 



This definition gives the correct mass of classical black hole solutions. The factor of 

12' 



(1 — Y2~e 2< ^) 2 does not affect this since it goes to 1 as x + — > oo but it is convenient for 
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what follows. For large negative x~ this expression for M(x~) gives the standard 
initial total energy of the infalling matter in gravitational collapse. As x~ increases 
the remaining mass will be reduced by an amount equal to the energy that has 
been radiated out in Hawking radiation up to that point. Using the equations of 
motion and the constraints one obtains 

fL = \ d+4>e -iPT™_ ) (3.18) 
ax A 

where T™_ is the - component of the matter energy momentum tensor with the 
contribution from the anomaly included. The factors in front of reflect the 
fact that the x + ,x~ coordinate system is not asymptotically Minkowskian. Thus 
(3.17) appears to be a sensible definition of the mass and it will clearly go to zero 
at the global event horizon if all the energy of the black hole is radiated away. 
M(x~) can be written as an integral as follows, 

oo 

MOO = Yxj dx+ ^ 1 ~^ h ~^ ' (3 - 19) 
o 

and by using the equations of motion and the constraints this can be worked into, 

OO 

M(x~) = J dx+ \j I - ^e 2 <t>e- 2p d-(l)T™+ . (3.20) 
o 

We will make use of this expression in section 5. 

We close this section by presenting some preliminary numerical results on black 
hole evolution in this 1+1 dimensional model. We consider a black hole formed by 
an incoming shock wave. When viewed on lines of constant x + the equations of 
motion (3.2) are ordinary differential equations for d + (fi and d + p, 

d-(d+<l>) ={1 - £r£*)d-(d +P ) , 

N N ^ 

2(1 - -e 2 *)d-{d + 4>) =(1 - -e 2 ^) (49^+0 + A 2 e 2 ") . 

This was used in reference [7] to obtain exact information about the solution just 
above the infall line and to establish the formation of the singularity at 

y|e 2 ^ = 1. The equations (3.21) can be solved exactly on the infall line because 
the functional coefficients in them are explicitly known there in terms of the linear 
dilaton background. But now the information about d+4> and d+p on the infall 
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line can be used to estimate the solution on the line x + = Xq + e + where e + is 
some small step size. Then the ordinary differential equations (3.21) can be solved 
numerically for d+(f) and d+p on that line, making it possible to take another step in 
the x + direction. Figure 3 shows typical results of evolving a gravitational collapse 
solution of the CGHS-equations in this way. The numerical integration cannot get 
through the singularity but since it lies on a space-like curve this does not affect 
the numerical solution in the physical region. The apparent horizon recedes as 
expected and it appears to approach the singularity asymptotically. 

We have also located some contours of constant (which are not shown in fig- 
ure 3) and evaluated the curvature scalar along them. The curvature does decrease 
along these contours as expected from (3.11). Numerical instability prevents us 
from following the contours very far into the x + direction and we can not at present 
determine from our data whether R goes to zero asymptotically. The hyperbolic 
nature of the CGHS-equations will make most numerical approaches unstable but 
one can certainly do better with a more sophisticated numerical method that the 
one we have employed here. 

4. Static Black Holes 

In this section we will discuss some properties of static black holes and speculate 
about their possible role as massive stable remnants of black hole evaporation. 
Some related issues are discussed in [12] We will present numerical evidence for 
the existence of static, non-radiating solutions of the CGHS-equations, which are 
candidates for stable 1+1 dimensional remnants. These solutions exist for any 
value of black hole mass and they have some interesting properties. In particular, 
the coupling is weak, Ng 2 « 1, everywhere in space if the mass is large, and thus 
it is possible that all further quantum corrections to these solutions remain small. 

At large positive x the static solutions rapidly tend to their classical counter- 
parts of the same mass but there is no energy in Hawking radiation in this region. 
For a massive black hole, M » NX, the classical geometry is well approximated 
until one approaches the position of the classical horizon, but there, in a metrically 
small region, the metric and dilaton change radically. In particular the curvature 
first increases to a maximum value, which depends on the mass, and then rapidly 
decreases, changes sign and goes to minus infinity at a finite proper distance, as 
one advances in the minus x direction. More surprisingly the dilaton reaches a 



★ These numerical solutions have also been obtained by J. Russo. While preparing this paper 
for publication we learned that B. Birnir, S. Giddings, J. Harvey and A. Strominger have 
found numerical static solutions to the CGHS-equations [13] and that they have also been 
studied by S. Hawking [14]. 
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maximum value at the same place as the curvature and thereafter the coupling 
strength decreases monotonically to zero as one approaches the curvature singu- 
larity. The maximum coupling is always smaller than the coupling in the classical 
solution at the horizon, i.e. g 2 < jj. This means that one might not expect the 
solution to receive any large corrections if M » NX. 

As mentioned in the introduction, massive remnants of this sort offer the pos- 
sibility of carrying off information about the quantum state of a black hole without 
leading to an infinite density of states at the Planck scale. The naked curvature 
singularity might pose problems of its own but these solutions are not ruled out 
on any a priori grounds. Later on we will consider the question of whether black 
hole evaporation can lead to one of these massive remnants as a final state. The 
answer turns out to be negative, at least in 1+1 dimensions, for we are able to 
show that the gravitational collapse solutions of the CGHS-equations cannot ap- 
proach massive static black holes asymptotically. This of course does not preclude 
the possibility of tunneling into one of these states at some stage in the evolution, 
but that process is presumably suppressed by large barrier factors for macroscopic 
black holes. 

Before we can rule out the static solutions as the final states, we need to explore 
some of their properties. We begin with a brief digression on the quantum theory of 
a free field in Rindler space [15,16]. This simple system shares important features 
with black holes both in four and two spacetime dimensions. 

The \x\ > \t\, x > wedge of flat two-dimensional Minkowski space may be 
described by Rindler coordinates R, 9 as follows, 

t =i?sinh^ , 

(4.1) 

x =i?cosh 9 , 
or x^ = ±Re ±0 . The Rindler metric is 



ds 2 = -R 2 d9 2 + dR 2 . (4.2) 

The Rindler Hamiltonian is conjugate to the time variable 9 and is simply the 
generator of Lorentz boosts in Minkowski space. At t — 0, 

oo 

H R = J dxT 00 (x)x. (4.3) 
o 

The Rindler Hamiltonian only acts on a Hilbert space defined in terms of fields 
living on the positive x-axis and a second Hilbert space is needed in order to 
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describe degrees of freedom on the negative x-axis. The Rindler Hamiltonian has a 
set of eigenvectors including a ground state. However, this ground state is not to be 
identified with the physical vacuum of Minkowski space, which has non-vanishing 
correlations between degrees of freedom on the positive and negative x-axis. These 
correlations are strongest at short distances and the Rindler ground state deviates 
violently from the physical vacuum very close to the origin. In fact, the expectation 
value of the Minkowski energy density diverges in the Rindler ground state at that 
point. Since the Minkowski vacuum is defined on the full x-axis it can not be a 
pure Rindler state but is rather a thermal density matrix, as is well known [15,16]. 

An analogous situation arises when a black hole is described in Kruskal co- 
ordinates. In that case the "Rindler ground state" is singular at the horizon. It 
describes a matter field configuration with divergent energy density which is ex- 
pected to react strongly back on the geometry at the horizon and therefore this state 
is usually not taken seriously. It nevertheless remains an open question whether 
there exists a corresponding zero temperature quantum state when gravitational 
corrections have been taken into account. It turns out that such a configuration 
can be found in our 1+1 dimensional model, as we shall shortly see, and we expect 
one to exist in 3+1 dimensions also. In fact, this is presumably the only static 
macroscopic field configuration in quantum gravity with a given mass. A finite 
temperature black hole in thermal equilibrium with its surroundings is not con- 
sistent because the thermal bath filling space would have infinite total energy. Of 
course, real black holes with finite Hawking temperature are not described by a 
static geometry but are formed in gravitational collapse. Hawking's original calcu- 
lation showed that gravitational collapse results in a black hole with a non-singular 
horizon and outgoing radiation with a thermal spectrum, but the black hole is not 
in thermal equilibrium with the rest of space. Whilst it is well established that 
black holes are formed with non-vanishing Hawking temperature and emit radia- 
tion it is very unclear how they subsequently evolve. It is an intriguing possibility 
that the back-reaction of Hawking radiation on the geometry could cause the evap- 
orating black hole to eventually settle into one of these massive zero temperature 
remnants, and we see no reason to reject this out of hand. We will, however, be 
able to show that this scenario is not realized in actual gravitational collapse, at 
least not in our 1+1 dimensional toy model. 

Let us first consider the energy momentum of quantum mechanical matter in a 
static black hole background ignoring back-reaction on the metric and dilaton. The 
static black hole solution (2.7) written in asymptotically Minkowskian coordinates 
r,C is 
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Far away from the black hole the metric is nearly flat and (p approaches the lin- 
ear dilaton background. In Euclidean space this solution is easily seen to be a 
coordinate transformation of Witten's semi-infinite cigar [17], the cigar tip being 
at C —> —oo. The Hawking effect can be described in terms of the conformal 
anomaly [4,5], and it is convenient to express it in terms of the Z field in (3.4). 
The curvature acts as a source term in the Z equation of motion, 

(4.5) 

The cigar is flat far away from the horizon and one finds a static Z field linear in 
the spatial coordinate, 

^/fc- (4-6) 

Inserting this into (3.5) gives a constant energy density, of Hawking radiation. 

The coefficient \J^r only involves the integrated curvature, or Euler number, of 
the cigar end, making it clear that Hawking temperature is independent of mass 
in 1 + 1 dimensions. There also exists a solution for Z which has vanishing energy 
density far away but that requires a point source of Z at the cigar tip and this 
causes the energy density to diverge at the horizon. This singular configuration is 
the zero temperature "Rindler ground state" . 

Neither the thermal black hole nor the zero temperature one remain consistent 
solutions when back- reaction is taken into account. As we mentioned before a 
thermal static configuration can not be consistent for it would have infinite total 
energy. On the other hand it is reasonable to look for a zero temperature static 
solution of the CGHS-equations with some finite mass (or ADM-energy). 

Static fields 0, p only depend on the spatial variable ( and satisfy the following 
set of ordinary differential equations in conformal gauge, 

0"=(1-|>)/, 

24 (4 7) 

(l-^)^=2(l-^)(0' 2 + A 2 e^), 



* Once we are willing to entertain such a singularity the temperature in the asymptotic region 
can of course have any value. 
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along with one constraint which can be written, 

0' 2 - P '4>' + ^P' 2 - \ 2 e 2p = . (4.8) 

The classical static black hole (4.4) is a solution of the equations without the Ne 2 ^ 
terms. Far away from the black hole the correction terms tend to zero and a given 
solution of the full equations matches onto a classical one with the same ADM- 
energy. Unfortunately we have not been able to find the general solution to the 
full set of non-linear equations (4.7) in closed form, but they can easily be handled 
numerically. In figures 5 to 6 we display some results of numerical integration. We 
plot both the value of e 2 ^ and the scalar curvature as a function of spatial position 
for af number of solutions of (4.7) with different black hole mass. Figure 4 graphs 
the coupling strength and curvature vs. position in the classical solution (4.4) for 
comparison. 

The figures clearly demonstrate the behavior alluded to at the beginning of this 
section. In particular, we see that the geometry and dilaton approach a classical 
black hole solution at large distances but before the classical horizon is reached 
the curvature reaches a maximum and the decreases to — oo. Simultaneously the 
coupling strength e 2< ^ reaches a maximum and then tends to zero as the curva- 
ture tends to — oo. The maximum coupling is always less than the critical value, 
N_ e 2ct>max <- i j-,^ approaches this value as M — > 0. The curvature reaches a maxi- 
mum positive value before it dips down to — oo. As M — > this maximum positive 
curvature diverges but at the same time the size of the region which deviates from 
the classical solution goes to zero. 

The solutions found above are the only static solutions of the CGHS-equations. 
As expected, no solutions exist which describe a static black hole in thermal equi- 
librium with its surroundings. 

A natural question to ask is whether a black hole formed in a physical process 
such as collapse can evolve into a massive zero temperature object described by 
these static solutions. We will now show that the CGHS-equations do not have 
solutions in which this occurs. That may not prevent quantum mechanical tunnel- 
ing into a massive zero temperature configuration at some stage in the evolution. 
However, we expect that to be heavily suppressed at least until the black hole has 
lost almost all its mass. 

We begin by noting that in the gravitational collapse of a shock wave the dilaton 
field increases monotonically on the infall line X — X Q . until it hits the singularity 
at ^e 2 ^ = 1. We will prove that this continues to hold on every line of constant 
x + , i.e. that it is not possible for the dilaton field to develop a local maximum 
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at j^e 2 ^ < 1 as is found in all M > static solutions. We know that on each 
such line e 2 ^ varies continuously from zero, as x~ — > — oo, to ^ at the singularity. 
(Recall that the singularity is a space-like curve which asymptotically approaches 
the light-like line, x~ = H, as shown in figure 2.) Thus if a local maximum occurs 
it must be accompanied by a local minimum closer to the singularity, as illustrated 
in figure 7. Call the points where the extrema occur x^(x + ) and At 
these points d-<fi = and from the second equation of motion in (3.2) it is seen 
that <9 + cL0 > there. It follows that the two curves x^(x + ) and x 2 {x + ) must 
approach one another as x + increases. However, since the (j) is monotonic on the 
infall line there must be a smallest value of x + where the extrema first occur and at 
which the equations of motion prevent x 1 and x 2 from separating 

it follows that <fi must remain monotonic. 

The local maxima of as a function of the spatial variable £ are also local 
maxima when viewed along constant x + lines. Therefore the CGHS-equations do 
not permit infalling matter to evolve to a massive remnant. The only static solution 
which has monotonic 4> is the M = solution. It seems highly plausible that the 
asymptotic endpoint of 1+1 dimensional black hole evolution is the M = state. 
Late in the evolution the system is expected to become highly quantum mechanical 
and the above argument does not rule out the possibility of forming light stable 
remnants. 

The approach to the vacuum state, as described by the CGHS-equations, is 
subtle and non-uniform. We believe that along any line of constant <p in the 
physical region, j^e 2 ^ < 1, the evaporating solution will eventually tend to the 
vacuum configuration. However, this cannot occur until the apparent horizon has 
receded beyond that value of <fi. As x + increases there will always be a small but 
diminishing region which remains inside the apparent horizon. In this region the 
solution is qualitatively different from the vacuum, for example constant contours 
are space-like there, but any given <fi contour will eventually emerge outside the 
horizon and become vacuum-like. It is unclear how much physical meaning should 
be attached to the Liouville region, y^e 2 ^ > 1. For our purposes the singularity 
can be taken as a boundary of spacetime. 
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5. The Fate of Quantum Information 



If stable remnants do not store the information that falls into a black hole then 
the only way to have a unitary quantum theory is to find a complete correlation 
between the state of the infalling matter which formed the black hole and the final 
state of the Hawking radiation. The fact that the final radiation is approximately 
thermal is in no way inconsistent with this. However, it is important to exhibit 
a mechanism for transferring the infalling information to the outgoing particles. 
Consider figure 8 which depicts an infalling system with some complexity forming 
a black hole and subsequently radiating away its energy. In order that the Hawking 
radiation carry away all the information it is necessary that the infalling matter 
has had all its information eliminated before it crosses the global event horizon at 
x~ = H. If the infalling matter has any memory of the initial state beyond this 
point then the Hawking radiation cannot be in a pure state since it would still be 
correlated to the infalling matter which has passed inside the event horizon. It is 
not sufficient for the infalling matter to imprint its state on the outgoing radiation 
like DNA imprinting itself on RNA. The original information in the "DNA" must 
also be destroyed if the "RNA" is to be in a pure state. At first sight this seems 
obviously impossible because the /-waves satisfy free field equations of motion and 
any structure that they have is maintained all along the infall line. In fact a free 
massless left-moving / field has an infinite number of conserved quantities which 
characterize the exact shape of the wave. 

On the other hand, since all the energy of the system is emitted as Hawking 
radiation the total energy of the infalling matter must be zero when it arrives at 
x~ = H. One might wonder whether all the apparent information carried by the 
/ fields is real. At the moment we have no definitive answer to this important 
question but we can exhibit a similar system in which the apparent information in 
an infalling free field is illusory and in which all real information is radiated in an 
analog of Hawking radiation. 

Consider the 1+1 dimensional Schwinger model* with a coupling constant 
which depends on position like 

9 2 = e 2x , (5.1) 

where x is the spatial coordinate. This can be thought of as electrodynamics in 
the linear dilaton vacuum. The Lagrangian is 

L = -j4r-^F^ + i^{d» - . (5.2) 



* 1+1 dimensional electrodynamics has been considered in a different context of black hole 
physics in reference [18]. 
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The analog for vector fields of the conformal gauge is the light-cone gauge A- = 0. 
In this gauge the Lagrangian becomes 

1 d d 

L = -^( d - A +) 2 + ^g—^L + iM-Qp: - i A +)i>R ■ (5-3) 

The field vpi describing "infalling" fermions appears to be a free field. A message 
could be encoded in sequence of particles and antiparticles of this type. Naively, 
any such information in ip^ is conserved and penetrates indefinitely into the strong 
coupling region. However, this is incorrect. To see this let us integrate out the i/j 
field. As is well known this contributes a one-loop correction to the effective action 
of F^, which is given by, 

l^J d x e^F, v ^F Xa . (5.4) 

Note the similarity with the Liouville action (3.1) obtained by integrating out the 
matter fields fi. This action is non-local but we can introduce a Z field as before. 
The effective Lagrangian then becomes 

L = \ 8 ^ z + 7S e ^ z ~ ■ (5 ' 5) 

Now integrate the gauge field to obtain 

L Z = \d,Zd»Z + 9 ^Z\ (5.6) 

This is the bosonized form of the Schwinger model in which the axial anomaly 
appears at the classical level. In terms of Z the charge current is given by e^ v d v Z. 

It is evident from (5.6) that information cannot penetrate indefinitely into the 
strong coupling region x — > — oo. The effective mass of the bosonizing field Z 
grows to infinity with the coupling so that any incoming Z wave of finite energy 
is completely reflected. This simply indicates the fact that electric current cannot 
flow into the infinitely strongly coupled region. 

The possibility of information penetrating the strong coupling region depends 
in detail on the model. For example if more than one fermion species carries charge 
then bosonization reveals one "massive" boson Z with mass ~ g 2 (x) and a set of 
massless bosons which can carry information toward x — > — oo. On the other hand 
if all currents are independently gauged then all the boson fields are reflected. 
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It is not clear whether or not a theory can be constructed in which information 
is always drained out of the / fields before they slip behind the global event horizon. 
It may be that an arbitrary theory, for example one which has global symmetries, 
cannot satisfy this. Perhaps it is only in some theories, or even only one, that 
black hole evolution can be made consistent with conservation of information. 

A possible candidate for such a theory in 1+1 dimensions is obtained if we 
completely replace the / fields in the CGHS-model by the field Z. This way 
the matter carries no information which is uncoupled to gravity. The question 
is whether in such a theory a configuration on the global event horizon which 
carries zero total energy must also carry no information. For example if T+ + in 
our expression for the energy (3.20) were a positive quantity as it would be in a 
classical theory then vanishing energy would require the Z field be trivial on the 
event horizon and therefore carry no information. Unfortunately energy density 
does not have to be positive in quantum theory. 

To begin to address this question consider an incoming configuration described 
asymptotically at large negative x~ by some left-moving positive energy density 
T ++ (x + ). It is always possible to find a Z field which carries the same energy 
distribution, 

T ++ (x+) = \v + ZV + Z + ^ V\Z . (5.7) 

Along a given line of fixed x~ one can find a coordinate transformation x + — > 
y + (x + ) which renders the metric flat along that line. The new coordinates are still 
conformal and (5.7) can be written as 

T ++ (*+) = _ e -v^(_) e ^ (5.8) 

along the line in question. If we write -0 = e^^ z then (5.8) takes a familiar 
form. It is the time independent Schrodinger equation for a zero energy state in a 
potential V(x + ) = ^T ++ (x + ). We can impose a boundary condition that Z = 
or tp — 1 in the linear dilaton vacuum at x + < Xq. If we assume that T ++ vanishes 
for large enough x + then the solution of the Schrodinger equation will be linear in 
that region. If the sign of the linear term is negative then ip will have to go through 
zero at some point, and Z becomes singular and ill-defined there. The sign depends 
on the potential V . For example, if J V(x + ) dx + < a zero energy solution of the 
Schrodinger equation always has a node somewhere. This means that an incoming 
asymptotic configuration with negative total energy would necessarily be singular 
and the only non-singular zero energy state has T ++ = 0. 
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In the case at hand it is not the incoming total energy which vanishes but 
rather the energy on the global event horizon. This condition does not appear 
to be strong enough to force Z to be vacuum-like along the entire global horizon. 
To see that introduce a y + coordinate system with flat metric on the global event 
horizon, x~ = H . The expression for the (vanishing) total energy (3.20) becomes 



/ 



dx+d-^1 - ^ e2 ^~ 1 (^+) 2 ^ = 0- ( 5 - 9 ) 



In general the weight factor d-<J)y 1 — y|e 2 ^ is not constant and that prevents us 
from concluding that Z must be trivial. In fact the weight factor increases with 
increasing x + and this might allow T++ to be positive on and near the infall line 
but to take a smaller compensating negative value at larger x + . Unless the back- 
reaction on the geometry provides more stringent constraints on T ++ on the global 
horizon then even in this "minimal" theory all information will not be drained out 
of the incoming Z field. 
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